KLRB1 defines an activated phenotype of CD4+ T cells and shows significant upregulation in patients with primary Sjögren's syndrome.

OBJECTIVE: This study aimed to investigate the role of KLRB1 (CD161) in human CD4+ T cells and elucidate its significance in primary Sjögren's syndrome (pSS). METHODS: Peripheral blood samples from 37 healthy controls and 44 pSS patients were collected. The publicly available single-cell RNA-Seq data from pSS patient PBMCs were utilized to analyse KLRB1 expression in T cells. KLRB1-expressing T lymphocyte subset proportions in pSS patients and healthy controls were determined by flow cytometry. CD25, Ki-67, cytokine secretion, and chemokine receptor expression in CD4+ KLRB1+ T cells were detected and compared with those in CD4+ KLRB1- T cells. Correlation analysis was conducted between KLRB1-related T-cell subsets and clinical indicators. ROC curves were generated to explore the diagnostic potential of KLRB1 for pSS. RESULTS: KLRB1 was significantly upregulated following T-cell activation, and Ki-67 and CD25 expression was significantly greater in CD4+ KLRB1+ T cells than in CD4+ KLRB1- T cells. KLRB1+ CD4+ T cells exhibited greater IL-17A, IL-21, IL-22, and IFN-γ secretion upon stimulation, and there were significantly greater proportions of CCR5+, CCR2+, CX3CR1+, CCR6+, and CXCR3+ cells among CD4+ KLRB1+ T cells than among CD4+ KLRB1- T cells. Compared with that in HCs, KLRB1 expression in CD4+ T cells was markedly elevated in pSS patients and significantly correlated with clinical disease indicators. CONCLUSION: KLRB1 is a characteristic molecule of the CD4+ T-cell activation phenotype. The increased expression of KLRB1 in the CD4+ T cells of pSS patients suggests its potential involvement in the pathogenesis of pSS and its utility as an auxiliary diagnostic marker for pSS.

FCRL3 expression is upregulated and closely correlates with TIGIT expression in regulatory T cells of patients with systemic lupus erythematosus.

Using data from single-cell RNA sequencing and flow cytometry, we initially examined the expression of FCRL3, finding it to be elevated and positively associated with TIGIT expression in the regulatory T cells of patients with systemic lupus erythematosus. This also suggests that the co-expression of FCRL3 and TIGIT warrants further attention.

Systematic characterization of p53-regulated long non-coding RNAs across human cancers reveals remarkable heterogeneity among different tumor types.

The p53 tumor suppressor protein, a sequence specific DNA binding transcription factor, regulates the expression of a large number of genes, in response to various forms of cellular stress. While the protein coding target genes of p53 have been well studied, less is known about its role in regulating long non-coding genes and their functional relevance to cancer. Here we report the genome-wide identification of a large set (>1000) of long non-coding RNAs(lncRNAs) that are putative p53 targets in a colon cancer cell line and in human patient datasets from five different common types of cancer. These lncRNAs have not been annotated by other studies of normal unstressed systems. In the colon cancer cell line a high proportion of these lncRNAs are uniquely induced by different chemotherapeutic agents that activate p53, while others are induced by more than one agent tested. Further, subsets of these lncRNAs independently predict overall and disease-free survival of patients across the five different common cancer types. Interestingly, both genetic alterations and patient survival associated with different lncRNAs are unique to each cancer tested, indicating extraordinary tissue-specific variability in the p53 non-coding response. The newly identified non-coding p53 target genes have allowed us to construct a classifier for tumor diagnosis and prognosis. Implications: Our results not only identify myriad p53-regulated lncRNAs, they also reveal marked drug-induced, as well as tissue- and tumor-specific heterogeneity in these putative p53 targets and our findings have enabled the construction of robust classifiers for diagnosis and prognosis.

Increased circulating regulatory T cells and decreased follicular T helper cells are associated with colorectal carcinogenesis.

Objective: Colorectal cancer (CRC) is the third most prevalent cancer worldwide and is associated with high morbidity and mortality rates. Colorectal carcinogenesis occurs via the conventional adenoma-to-carcinoma and serrated pathways. Conventional T helper (Th) and innate lymphoid cells (ILCs) play vital roles in maintaining intestinal homeostasis. However, the contribution of these two major lymphoid cell populations and their associated cytokines to CRC development is unclear. Therefore, we aimed to analyze peripheral lymphocyte profiles during colorectal carcinogenesis. Methods: We collected 86 blood samples concurrently, and pathologists confirmed the presence of various pathological conditions (i.e., HPs, adenoma, and carcinoma) using hematoxylin and eosin staining. Ten healthy donors were recruited as healthy controls (HCs) from the physical examination center. We performed flow cytometry on peripheral blood mononuclear cells collected from patients with various pathological conditions and the HCs, and cytokines (interleukin-2, interleukin-4, interleukin-5, interleukin-13, interleukin-17A, interleukin-17F, interleukin-22, interferon-γ, and tumor necrosis factor-α) were quantified. We also analyzed the published single-cell RNA sequence data derived from tissue samples from different stages of colorectal carcinogenesis. Results: The cytokine response in peripheral CD4+ T cells was upregulated during the carcinoma process. The frequency of peripheral regulatory T cells (Tregs) increased in the adenoma and carcinoma stages. While the T follicular helper (Tfh) cell proportion was downregulated in the adenoma and carcinoma processes. Thus, Th cell subsets, especially Tregs and Tfh cells, were involved in colonic diseases. Moreover, the immunological profile characteristics in the HPs were clarified. Conclusion: We comprehensively analyzed circulating ILCs and adaptive T-cell lymphocyte subtypes in colorectal carcinoma progression. Our results show the immunological profile characteristics and support the involvement of Th subsets, especially Treg and Tfh cell populations, in colonic diseases. These findings significantly enhance our understanding of the immune mechanisms underlying CRC and its precancerous lesions. Further investigation of the Treg and Tfh cells' function in colorectal disease development will provide potential therapeutic targets for monitoring and preventing CRC development.

Identification of rat Vstm1 with conservative anti-inflammatory effect between rat and human homologs.

Human VSTM1 (also known as SIRL1) is an inhibitory immune checkpoint receptor involved in leukocyte activation. Identification of the homologous genes in other species, such as mice and rats, will undoubtedly contribute to functional studies and clinical applications. Here, we successfully cloned the Vstm1 gene in rats, as supported by high-throughput sequencing data. However, Vstm1 is degenerated to a pseudogene in the mouse genome. Rat Vstm1 mRNA contains a complete open reading frame (ORF) of 630 nucleotides encoding 209 amino acids. Rat Vstm1 is highly expressed in bone marrow, especially in granulocytes. The expression levels of Vstm1 gradually increase with the development of granulocytes in bone marrow but are downregulated in response to inflammatory stimuli. Rat VSTM1 does not have an immunoreceptor tyrosine-based inhibitory motif (ITIM), however, it shows a conservative function of inflammatory inhibition with human VSTM1, and both are anti-correlated with many inflammatory cytokines, such as IL-1α and TNF-α. In bone marrow-derived macrophages (BMDMs), either rat or human VSTM1 suppressed the secretion of inflammatory cytokines in response to LPS stimulation. Further analysis in lung cancer microenvironment revealed that VSTM1 is mainly expressed in myeloid cells, anti-correlated with inflammatory cytokines and associated with tumor development and metastasis.

Exosomal secreted SCIMP regulates communication between macrophages and neutrophils in pneumonia.

In pneumonia, the deficient or delayed pathogen clearance can lead to pathogen proliferation and subsequent overactive immune responses, inducing acute lung injury (ALI). While screening human genome coding genes using our peripheral blood cell chemotactic platform, we unexpectedly find SLP adaptor and CSK interacting membrane protein (SCIMP), a protein with neutrophil chemotactic activity secreted during ALI. However, the specific role of SCIMP in ALI remains unclear. In this study, we investigate the secretion of SCIMP in exosomes (SCIMPexo) by macrophages after bacterial stimulation, both in vitro and in vivo. We observe a significant increase in the levels of SCIMPexo in bronchoalveolar lavage fluid and serum of pneumonia patients. We also find that bronchial perfusion with SCIMPexo or SCIMP N-terminal peptides increases the survival rate of the ALI model. This occurs due to the chemoattraction and activation of peripheral neutrophils dependent on formyl peptide receptor 1/2 (FPR1/2). Conversely, exosome suppressors and FPR1/2 antagonists decrease the survival rate in the lethal ALI model. Scimp-deficient and Fpr1/2-deficient mice also have lower survival rates and shorter survival times than wild-type mice. However, bronchial perfusion of SCIMP rescues Scimp-deficient mice but not Fpr1/2-deficient mice. Collectively, our findings suggest that the macrophage-SCIMP-FPRs-neutrophil axis plays a vital role in the innate immune process underlying ALI.

Upregulation of CX3CR1 expression in circulating T cells of systemic lupus erythematosus patients as a reflection of autoimmune status through characterization of cytotoxic capacity.

OBJECTIVE: This study investigated CX3CR1 expression in human peripheral blood T lymphocytes and their subsets, exploring changes in SLE patients and its diagnostic potential. METHODS: Peripheral blood samples from 31 healthy controls and 50 SLE patients were collected. RNA-Seq data from SLE patient PBMCs were used to analyze CX3CR1 expression in T cells. Flow cytometry determined CX3CR1-expressing T lymphocyte subset proportions in SLE patients and healthy controls. Subset composition and presence of GZMB, GPR56, and perforin in CX3CR1+ T lymphocytes were analyzed. T cell-clinical indicator correlations were assessed. ROC curves explored CX3CR1's diagnostic potential for SLE. RESULTS: CX3CR1+CD8+ T cells exhibited higher GPR56, perforin, and GZMB expression than other T cell subsets. The proportion of CX3CR1+ was higher in TEMRA and lower in Tn and TCM. PMA activation reduced CX3CR1+ T cell proportions. Both RNA-Seq and flow cytometry revealed elevated CX3CR1+ T cell proportions in SLE patients. Significantly lower perforin+ and GPR56+ proportions were observed in CX3CR1+CD8+ T cells in SLE patients. CX3CR1+ T cells correlated with clinical indicators. CONCLUSION: CX3CR1+ T cells display cytotoxic features, with heightened expression in CD8+ T cells, particularly in adult SLE patients. Increased CX3CR1 expression in SLE patient T cells suggests its potential as an adjunctive diagnostic marker for SLE.

CMTM6 promotes hepatocellular carcinoma progression through stabilizing ß-catenin.

CMTM6, a regulator of PD-L1 stability, has been implicated in the development of various cancers. However, the expression and role of CMTM6 in hepatocellular carcinoma (HCC) remains controversial. Our study revealed a negative correlation between CMTM6 expression and HCC prognosis through bioinformatics analysis and immunofluorescence staining. CMTM6 expression was also positively associated with alpha-fetoprotein (AFP) levels, supporting its potential as a prognostic marker for HCC. Using Cmtm6 knockout mice, we found that Cmtm6 deficiency inhibited HCC formation and cell proliferation in primary liver cancer models induced by DEN and DEN/CCl4. In HCC cell lines, CMTM6 promoted cell proliferation and interacted with ß-catenin, stabilizing it by preventing ubiquitination. In conclusion, our study suggested that CMTM6 upregulation promotes HCC cell proliferation through the ß-catenin pathway, making it a potential therapeutic target for HCC treatment.

Omic horizon expression: a database of gene expression based on RNA sequencing data.

BACKGROUND: Gene expression profiles have important significance for gene expression characteristics and further functional studies. More attention has been given to the expression databases in humans and mice, but less attention has been given to rats, while rat models also play an irreplaceable role in biomedical experiments. RESULTS: To depict the rat gene expression profiles in mRNA expression levels, we analyzed over 2,700 RNA sequencing (RNA-Seq) samples from 48 tissues, 40 primary cell types and 25 cell lines; and then mapped them to the latest version of the rat genome reference, mRatBN7.2. Based on these datasets and reanalysis, we constructed a new database, the Omic Horizon Expression Database ( http://immudb.bjmu.edu.cn/expression.html ), which allows expressional profile query of over 25,000 rat genes based on non-redundant gene symbols. The database supports requests using gene symbols (or alias), Ensemble and Entrez gene IDs. Gene expression profiles can be queried in three categories: tissues, primary cells and cell lines. Application examples including expression profiling and comparison, as well as identification of novel rat genes, were illustrated to show the utility of the database. CONCLUSIONS: As an omic resource, the Omic Horizon Expression Database provides horizons of gene expression profiles across various tissues and cells, which greatly facilitates the identification of rat genes as well as functional clues.

A pan-cancer analysis implicates human NKIRAS1 as a tumor-suppressor gene.

The NF-κB family of transcription factors and the Ras family of small GTPases are important mediators of proproliferative signaling that drives tumorigenesis and carcinogenesis. The κB-Ras proteins were previously shown to inhibit both NF-κB and Ras activation through independent mechanisms, implicating them as tumor suppressors with potentially broad relevance to human cancers. In this study, we have used two mouse models to establish the relevance of the κB-Ras proteins for tumorigenesis. Additionally, we have utilized a pan-cancer bioinformatics analysis to explore the role of the κB-Ras proteins in human cancers. Surprisingly, we find that the genes encoding κB-Ras 1 (NKIRAS1) and κB-Ras 2 (NKIRAS2) are rarely down-regulated in tumor samples with oncogenic Ras mutations. Reduced expression of human NKIRAS1 alone is associated with worse prognosis in at least four cancer types and linked to a network of genes implicated in tumorigenesis. Our findings provide direct evidence that loss of NKIRAS1 in human tumors that do not carry oncogenic RAS mutations is associated with worse clinical outcomes.

Loss of the Novel Myelin Protein CMTM5 in Multiple Sclerosis Lesions and Its Involvement in Oligodendroglial Stress Responses.

This study comprehensively addresses the involvement of the protein CKLF-like Marvel transmembrane domain-containing family member 5 (CMTM5) in the context of demyelination and cytodegenerative autoimmune diseases, particularly multiple Sclerosis (MS). An observed reduction in CMTM5 expression in post-mortem MS lesions prompted further investigations in both in vitro and in vivo animal models. In the cuprizone animal model, we detected a decrease in CMTM5 expression in oligodendrocytes that is absent in other members of the CMTM protein family. Our findings also confirm these results in the experimental autoimmune encephalomyelitis (EAE) model with decreased CMTM5 expression in both cerebellum and spinal cord white matter. We also examined the effects of a Cmtm5 knockdown in vitro in the oligodendroglial Oli-neu mouse cell line using the CRISPR interference technique. Interestingly, we found no effects on cell response to thapsigargin-induced endoplasmic reticulum (ER) stress as determined by Atf4 activity, an indicator of cellular stress responses. Overall, these results substantiate previous findings suggesting that CMTM5, rather than contributing to myelin biogenesis, is involved in maintaining axonal integrity. Our study further demonstrates that the knockdown of Cmtm5 in vitro does not modulate oligodendroglial responses to ER stress. These results warrant further investigation into the functional role of CMTM5 during axonal degeneration in the context of demyelinating conditions.

Evaluating imputation methods for single-cell RNA-seq data.

BACKGROUND: Single-cell RNA sequencing (scRNA-seq) enables the high-throughput profiling of gene expression at the single-cell level. However, overwhelming dropouts within data may obscure meaningful biological signals. Various imputation methods have recently been developed to address this problem. Therefore, it is important to perform a systematic evaluation of different imputation algorithms. RESULTS: In this study, we evaluated 11 of the most recent imputation methods on 12 real biological datasets from immunological studies and 4 simulated datasets. The performance of these methods was compared, based on numerical recovery, cell clustering and marker gene analysis. Most of the methods brought some benefits on numerical recovery. To some extent, the performance of imputation methods varied among protocols. In the cell clustering analysis, no method performed consistently well across all datasets. Some methods performed poorly on real datasets but excellent on simulated datasets. Surprisingly and importantly, some methods had a negative effect on cell clustering. In marker gene analysis, some methods identified potentially novel cell subsets. However, not all of the marker genes were successfully imputed in gene expression, suggesting that imputation challenges remain. CONCLUSIONS: In summary, different imputation methods showed different effects on different datasets, suggesting that imputation may have dataset specificity. Our study reveals the benefits and limitations of various imputation methods and provides a data-driven guidance for scRNA-seq data analysis.

Corrigendum: A theoretical framework of immune cell phenotypic classification and discovery.

[This corrects the article DOI: 10.3389/fimmu.2023.1128423.].

G Protein-Coupled Receptor 56 Characterizes CTLs and Reflects the Progression of Lung Cancer Patients.

CTLs play important roles in host immune responses to tumors. CD4 CTLs are characterized by their ability to secrete cytotoxic effector molecules, such as granzyme B and perforin, and kill target cells in a MHC class II-restricted manner. However, the cell surface markers of CD4 CTLs remain unknown, which hinders their separation and research on their function. In this study, we performed a bioinformatics analysis and experimental validation that revealed that G protein-coupled receptor 56 (GPR56) is a cell surface marker that can be used to characterize CD4 CTLs. We found that GPR56 and granzyme B were coexpressed in extremely high levels in human peripheral blood T cells, and that anti-GPR56 stimulation significantly upregulated the expression of granzyme B in both CD4+GPR56+ and CD8+GPR56+ T cells. These findings suggest that GPR56 expression and the GPR56 signaling pathway could contribute directly to the toxic function of either CD4+ or CD8+ T cells. We also used GPR56 as a biomarker to investigate the clinical significance of CD4 CTLs. GPR56+ T cell levels were increased in patients with lung cancer, and GPR56 expression was significantly correlated with lung cancer progression. A further analysis revealed an increase in exhausted cell states in lung cancer patients because of upregulation of programmed cell death protein 1 expression in GPR56+ T cells. The findings of this study suggest that GPR56 characterizes the cytotoxic states of either CD4+ or CD8+ T cells.

B-Cell Receptor Features and Database Establishment in Recovered COVID-19 Patients by Combining 5'-RACE with PacBio Sequencing.

BACKGROUND: Antibodies induced by viral infection can not only prevent subsequent virus infection, but can also mediate pathological injury following infection. Therefore, understanding the B-cell receptor (BCR) repertoire of either specific neutralizing or pathological antibodies from patients convalescing from Coronavirus disease 2019 (COVID-19) infection is of benefit for the preparation of therapeutic or preventive antibodies, and may provide insight into the mechanisms of COVID-19 pathological injury. METHODS: In this study, we used a molecular approach of combining 5' Rapid Amplification of cDNA Ends (5'-RACE) with PacBio sequencing to analyze the BCR repertoire of all 5 IgH and 2 IgL genes in B-cells harvested from 35 convalescent patients after severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. RESULTS: We observed numerous BCR clonotypes within most COVID-19 patients, but not in healthy controls, which validates the association of the disease with a prototypical immune response. In addition, many clonotypes were found to be frequently shared between different patients or different classes of antibodies. CONCLUSIONS: These convergent clonotypes provide a resource to identify potential therapeutic/prophylactic antibodies, or identify antibodies associated with pathological effects following infection with SARS-CoV-2.

Mapping Resident Immune Cells in the Murine Ocular Surface and Lacrimal Gland by Flow Cytometry.

BACKGROUND: The ocular surface and lacrimal gland have a frontline position in mucosal immunology. However, there have been few updates to the immune cell atlas of these tissues in recent years. PURPOSE: To map the immune cells in murine ocular surface tissues and lacrimal gland. METHODS: Central and peripheral corneas, conjunctiva, and lacrimal gland were dissociated into single cell suspensions, followed by flow cytometry. Discrepancy of immune cells between the central and peripheral corneas was compared. In the conjunctiva and lacrimal gland, myeloid cells were clustered by tSNE and FlowSOM based on the expression of F4/80, Ly6C, Ly6G, and MHC II. ILCs, type 1 immune cells, and type 3 immune cells were analyzed. RESULTS: The number of immune cells in peripheral corneas was about 16 folds of that in central corneas. B cells accounted for 8.74% of immune cells in murine peripheral corneas. In the conjunctiva and lacrimal gland, most myeloid cells tended out to be monocytes, macrophages, and classical dendritic cells (cDCs). ILC3 were 6.28% and 3.63% of ILCs in the conjunctiva and lacrimal gland, respectively. Th1, Tc1, and NK cells were predominant type 1 immune cells. γδ T17 cells and ILC3 outnumbered Th17 cells among type 3 T cells. CONCLUSION: B cells resident in murine corneas were reported for the first time. Additionally, we proposed a strategy of clustering myeloid cells to better understand their heterogeneity in the conjunctiva and lacrimal gland based on tSNE and FlowSOM. Furthermore, we identified the ILC3 in the conjunctiva and lacrimal gland for the first time. Compositions of type 1 and type 3 immune cells were summarized. Our study provides a fundamental reference and novel insights for ocular surface immune homeostasis and diseases.

A theoretical framework of immune cell phenotypic classification and discovery.

Immune cells are highly heterogeneous and show diverse phenotypes, but the underlying mechanism remains to be elucidated. In this study, we proposed a theoretical framework for immune cell phenotypic classification based on gene plasticity, which herein refers to expressional change or variability in response to conditions. The system contains two core points. One is that the functional subsets of immune cells can be further divided into subdivisions based on their highly plastic genes, and the other is that loss of phenotype accompanies gain of phenotype during phenotypic conversion. The first point suggests phenotypic stratification or layerability according to gene plasticity, while the second point reveals expressional compatibility and mutual exclusion during the change in gene plasticity states. Abundant transcriptome data analysis in this study from both microarray and RNA sequencing in human CD4 and CD8 single-positive T cells, B cells, natural killer cells and monocytes supports the logical rationality and generality, as well as expansibility, across immune cells. A collection of thousands of known immunophenotypes reported in the literature further supports that highly plastic genes play an important role in maintaining immune cell phenotypes and reveals that the current classification model is compatible with the traditionally defined functional subsets. The system provides a new perspective to understand the characteristics of dynamic, diversified immune cell phenotypes and intrinsic regulation in the immune system. Moreover, the current substantial results based on plasticitomics analysis of bulk and single-cell sequencing data provide a useful resource for big-data-driven experimental studies and knowledge discoveries.

Role of FABP5 in T Cell Lipid Metabolism and Function in the Tumor Microenvironment.

To evade immune surveillance, tumors develop a hostile microenvironment that inhibits anti-tumor immunity. Recent immunotherapy breakthroughs that target the reinvigoration of tumor-infiltrating T lymphocytes (TIL) have led to unprecedented success in treating some cancers that are resistant to conventional therapy, suggesting that T cells play a pivotal role in anti-tumor immunity. In the hostile tumor microenvironment (TME), activated T cells are known to mainly rely on aerobic glycolysis to facilitate their proliferation and anti-tumor function. However, TILs usually exhibit an exhausted phenotype and impaired anti-tumor activity due to the limited availability of key nutrients (e.g., glucose) in the TME. Given that different T cell subsets have unique metabolic pathways which determine their effector function, this review introduces our current understanding of T cell development, activation signals and metabolic pathways. Moreover, emerging evidence suggests that fatty acid binding protein 5 (FABP5) expression in T cells regulates T cell lipid metabolism and function. We highlight how FABP5 regulates fatty acid uptake and oxidation, thus shaping the survival and function of different T cell subsets in the TME.

A new strategy: identification of specific antibodies for neutralizing epitope on SARS-CoV-2 S protein by LC-MS/MS combined with immune repertoire.

Although the SARS-CoV-2 vaccine has been widely used worldwide, not all individuals can produce neutralization antibodies, so it is still urgent to find and prepare neutralization antibodies for COVID-19 prevention or treatment. In this study, we created a new strategy to effectively obtain neutralizing antibodies or complementary determining region 3 (CDR3) of neutralizing antibodies against SARS-CoV-2. We first predicted and synthesized several B cell epitopes on RBD and adjacent RBD of S protein, then the B cell epitopes were used to prepare affinity chromatography columns respectively and purify the binding IgG from serum samples of convalescent COVID-19 patients. After these IgGs were identified to have neutralizing activity, the peptide sequences of the antigen-binding regions (variable region) of neutralizing antibodies were analyzed by protein mass spectrometry. Subsequently, the B cells from the same individual were sorted and used to obtain their full BCR repertoire by 5' RACE combined with high-throughput of PacBio sequencing method. Then, the peptide sequence of neutralizing antibody variable region by protein mass spectrometry was mapped to the full BCR repertoire and found the full variable region sequence of neutralizing antibodies. Finally, we obtained and synthesized numerous CDR3 peptides of neutralizing antibodies to confirm the neutralizing activity for SARS-CoV-2 infection. Our results indicate that the novel scheme will be suitable for rapid screening of neutralizing antibodies, including screening neutralizing antibodies against SARS-CoV-2 and other pathogenic microorganisms.

ImmuMethy, a database of DNA methylation plasticity at a single cytosine resolution in human blood and immune cells.

Differential DNA methylation is a feature of numerous physiological and pathological processes. However, the extent to which single-base cytosine methylation modifies cellular responses to various stimuli has not been well characterized. In this study, we carried out a systematic analysis of methylome data derived from human blood and immune cells and constructed the ImmuMethy database. ImmuMethy allows interrogation of DNA methylation plasticity (MPL) at the single cytosine level. MPL, which refers to the variability of DNA methylation, is quantitatively measured in multiple ways, such as quartiles and standard deviations. ImmuMethy comprises over 36 000 samples from the Human Methylation450 and MethylationEPIC BeadChips platforms and provides multiple applications, such as an overview of methylation status and plasticity, differential methylation analysis, identification of methylation markers and sample stratification. An analysis of all datasets revealed that DNA methylation is generally stable, with minimal changes in beta values. This further supports the characteristics of DNA methylation homeostasis. Based on the beta value distribution, we identified three types of methylation sites: methylation tendency sites, unmethylation tendency sites and dual tendency or nonbiased methylation sites. These sites represent different methylation tendentiousness of DNA methylation across samples. The occurrence of multiple methylation tendencies in a site means split methylation, which generally corresponds to high MPL. Inverted methylation tendencies from methylation tendency sites to unmethylation tendency sites, or vice versa, represent strong differential methylation in response to conditions. All these sites can be identified in ImmuMethy, making it a useful tool for omics-based data-driven knowledge discovery. Database URL: http://immudb.bjmu.edu.cn/immumethy/.